Non-Equilibrium Dynamics with Fermions on a Lattice in Space and Time* 
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We consider the dynamics of the 1 + 1 dimensional abelian Higgs model with axially coupled 
fermions, in the large Nf limit, on a lattice in space and real-time. We allow for inhomogeneous 
classical Bose fields. In order to deal with the lattice doublers, we use Wilson's lattice fermions. 
The lattice formulation leads to a stable integration algorithm. We illustrate the implementation 
with numerical results. 



I. INTRODUCTION 

Real-time dynamics of quantum fields plays an impor- 
tant role in the early universe, in heavy ion collisions 
and in condensed matter. For example, a detailed un- 
derstanding of the rate of sphaleron transitions at high 
temperature is relevant for baryogenesis. This is a non- 
perturbative problem, which should be treated preferably 
directly in real-time. 

However, as is usually the case, in order to do ac- 
tual calculations, approximations have to be made. The 
classical approximation is based on the observa- 

tion that at sufficiently high temperature the relevant 
long distance modes in the quantum theory behave clas- 
sically. Hence it makes sense to study the classical equa- 
tions of motion. The advantage is that it is (almost) 
straightforward to implement it numerically, and that 
e.g. sphaleron transitions are easy to find. However, it 
has become clear that Hard Thermal Loop (HTL) effects 
play an important role, and that it is necessary to incor- 
porate them in a classical field and/or classical particle 
approach 

Other possibilities are large N or Hartree approxima- 
tions [p|-pd|, which have been applied to e.g. (inflation- 
ary) scalar models and gauge theories. In this formula- 
tion the fields are split in a classical part and a quantum 
field, and they are coupled via semiclassical equations of 
motion, which are exact in the limit N — > oo. The advan- 
tage is that quantum effects, including the backreaction 
on the classical fields, are taken into account. However, in 
actual numerical calculations, the emphasis has been on 
the case of homogeneous classical fields . This ex- 

cludes inhomogeneous classical field configurations, such 
as the sphaleron. 

Since we want to improve both the classical approxi- 
mation by including quantum effects (such as HTL's in 
3 + 1 dimensions) and the large N approach by relaxing 
the classical field being homogeneous, we consider here 
the second approach with inhomogeneous classical fields. 



After a discussion of the model, we describe the effec- 
tive equations of motion. They are defined on a lattice in 
space and time, to have gauge invariance and numerical 
stability. We illustrate this with some numerical results. 



II. ABELIAN HIGGS MODEL WITH AXIALLY 
COUPLED FERMIONS 



We consider the abelian Higgs model in 1 + 1 dimen- 
sions, with axially coupled fermions. The (continuum) 
action is 
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It has a number of properties in common with the 
Standard Model. The bosonic vacuum is non-trivial, 
and is labeled with the Chern-Simons number C = 
— f dx 1 Ai / 2ir , being integer at a vacuum. The vacua 
are separated by a finite energy sphaleron barrier, where 
C is half integer. There is a global vector symmetry, 
■ijj — > e^tp, which leads to a classically conserved vec- 
tor current j*i — itp^/^ip. However, the corresponding 
charge, the fermion number Qf — J dx 1 ^, is not con- 
served in the quantum theory, because of the anomaly: 
Qf(t)-Q f (0) = ~{C{t) - C(0)). The bosonic part of 
this model has been used before to test the classical ap- 
proximation ^,|). 

Being in 1 + 1 dimensions has a few consequenses. The 
typical HTL expressions are sub dominant, so that these 
play a less important role j|. Also the divergences are 
less severe, only the Higgs self energy is divergent, and 
only due to Higgs and fermion one loop contributions. 
And it is much easier to deal with numerically than 3 + 1 
dimensional models. 
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III. ON THE LATTICE IN REAL-TIME 

In order to formulate the model on a lattice, it is con- 
venient to first put the model in a more suitable form. 
By performing charge conjugation on the right-handed 
fermions only, ipn — > (-0^C) T , ipL —> ipL, where C is the 
charge conjugation matrix, we write the model as one 
with a vector gauge symmetry. The fermionic part of 
the action becomes the usual Schwinger model with a 
Majorana-like Yukawa coupling 
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The anomalous current is now axial, jg = 1^7^75^, and 
the anomaly equation reads 

Q 5 (t) - Q 5 (0) = C(t) - C(0). 

Since the Yukawa term has become Majorana-like, it is 
convenient to work with a real 4-component Majorana 
field 

It is now straightforward to put this theory on a space- 
time lattice with lattice spacing <zo,ai = a, using Wil- 
son's lattice fermions, known from euclidean lattice field 
theory. Lattice equations of motion and currents associ- 
ated to global symmetries follow in the usual way. 



IV. EFFECTIVE EQUATIONS OF MOTION 

Effective equations of motion can be derived by consid- 
ering Nf fermion fields with Nf — > 00. After integrating 
out the fermions and rescaling <f> and v with U Nf, e and 
G with 1/ \fNf, and A with l/Nf, the resulting bosonic 
effective action is proportional to Nf, and the remain- 
ing bosonic path integral can be approximated with a 
saddle point expansion. In leading order this results in 
the semiclassical Maxwell-like equations of motion || . In 
this approximation, the Bose fields are treated as classi- 
cal fields, and the fermions, that are fully quantized, live 
in their background. The initial conditions for the origi- 
nal (unsealed) <f> and eA^ have evidently to be such that 
they are of order y/Nf. 

When the classical field is constrained to be homoge- 
neous, the quantum fluctuations are treated in practise 
with a mode function expansion |Io|JTl| l . Since we want to 
allow also inhomogencous Bose fields, we use a 'general- 
ized' mode function expansion, i.e. at t = 0, the fermion 
field is expanded in plane waves in order to be able to 
specify the initial quantum state, but for t > the space 
dependence is determined by the field equations. The re- 
sult is that the mode functions have both a momentum 
label p corresponding to the initial state, and a general 
space (and time) dependence, written as x = (t, x 1 ). 

We expand the field as 



)w p (x; a) + b P {a)w*{x; a)] 



with the initial condition 

w p (t = 0,x 1 ;a) = e ipxl w p (a). 

Wp{pi) is an eigenspinor of the Dirac hamiltonian at t = 0, 
and a = {+, — } labels the independent eigenspinors. Ex- 
pectation values of the time independent creation and an- 
nihilation operators bp^ (a) specify the quantum state at 
t = 0. We take either a vacuum or a thermal initial state, 
summarized by c p (a) = ([b p (a) , b p (a)]) = 1 — 2f p (a), 
with f p the Fermi-Dirac distribution. 

The bosonic equations are (in the temporal gauge 
A = 0) 

d'd A lx = e 2 (jl + (j}J), 
^ A lx = -e 2 0l+0X>), 
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The second equation is Gauss' law. The (primed) deriva- 
tives denote forward (backward) lattice derivatives, is 
the classical Higgs contribution to the current. 

The fermion current (jj? ) and force due to the fermions 
on the scalar field (F) can be written in terms of mode 
functions. The mode functions themselves obey the lat- 
tice Dirac equation in the presence of the Bose fields. 
This gives a closed set of equations. 

Also observables such as the fermion vector and axial 
charge are expressed in terms of the mode functions. As 
an example, we show here the axial charge density (j® ) 
on the lattice 
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Here p 2 is a matrix appearing in the Majorana descrip- 
tion. 

The equations of motion still contain the divergence 
of the quantum theory. Since we only integrated out 
the fermions, and we treat the scalar field completely 
classical, only the fermion loop divergence enters in the 
equations. (F x ) is a logarithmically divergent sum. It's 
divergence is canceled with the appropriate bare Vg. 



V. NUMERICAL RESULTS 

We have used the following initial conditions. The ini- 
tial fermion mode functions are eigenspinors of the Dirac 
hamiltonian in the presence of a vacuum configuration 
of Bose fields, i.e. A\ = 0, <j> = vr/\/2, and we take a 
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vacuum initial state for the creation and annihilation op- 
erators, i.e. a p (a) — 1. The Higgs field is initially homo- 
geneous, <f> — vr/ \/2, but for its time derivative, do<j) = ?r> 
only the low momentum modes have a non-zero ampli- 
tude. The high momentum modes of 7r are initially zero. 
This makes the classical field inhomogeneous. Both the 
Chern-Simon number and its time derivative have a non- 
zero value at t = 0. The initial configuration is chosen 
in such a way that it obeys Gauss' law. This completely 
specifies the spatial dependence of doAi. 
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FIG. 1. Qh, Qf and Q to t — Qf + Qh versus et. 
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FIG. 2. Total conserved charge Qtot versus et. 

The quantum divergence in the equations is treated in 
the following way. The renormalized Higgs expectation 
value vr is fixed, and then the corresponding bare vb 
is determined using the equation for the scalar field at 
t = 0. We checked that this leads to convergent physical 
results when the lattice spacing a is made smaller. 

Every dimensionful parameter is written in units of e. 
The parameters are A/e 2 = G/e — 0.25 and the renor- 
malized Higgs expectation value is v 2 R — 10. The physical 
size is eL = 3.2 (L — aN). Furthermore, the temporal 
lattice spacing ao/a = 0.01, and the number of lattice 



Gauss' law de- 
Q h {t) +Q f {t) 



points N = 64. 

We now discuss the numerical results, 
mands that the total charge Qtot(i) = 
is zero. Our initial conditions are such that both the 
Higgs charge, Qh(0), and the fermion charge, Q/(0), are 
separately zero. Under real-time evolution, having an in- 
homogeneous system, the individual charges do not stay 
zero, but the sum vanishes (up to machine precision). 
This is demonstrated in Figs. [l], @. Since we solve a large 
set of partial differential equations (in particular there is 
an equation for every mode function, and there are 2N 
mode functions), this is non-trivial. 
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FIG. 3. Bosonic energy, Eb/e, fermionic energy, Ef/e, and 
total energy E to t/e = (Eb + Ef)/e versus et. 
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FIG. 4. Total conserved energy Etot/e versus et. 

Also the total energy E to t = Eb + Ef is conserved, but 
not to such high accuracy as the charge. This is shown 
in Figs. |3| |i| The fermions are initialized in a vacuum 
state (Ef(0) = 0). All the initial energy is contained 
in the low momentum modes of the Bose fields. This is 
clearly a non-equilibrium situation, and energy is trans- 
ferred into the fermion degrees of freedom during time 
evolution. 
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Interesting observables are the Chern-Simons number 
C and the anomalous axial charge Q$. According to the 
anomaly equation, these should be equal during time evo- 
lution. C is a 'simple' observable in the sense that it is 
just the sum of A\ over all lattice points. On the contrary, 
is obtained by taking the sum over inner products be- 
tween all modes functions, see (0). Results are given m 
Fig. |[ It shows that Q5 follows C at early times, in 
agreement with the anomaly equation. For later times 
they start to deviate. This is an artefact of the lattice 
discretization, which is reduced by taking smaller lattice 
spacing. In Fig. @ we show the Higgs order parameter. 



in space and time. Because of the presence of 2N spinor 
mode functions that are space dependent, the method is, 
especially in 3 + 1 dimensions, numerically demanding. 

The 1 + 1 dimensional model is an excellent toy 
model to test the possibility of numerical calculations in 
a non-homogeneous, non-equilibrium context, including 
anomalous fermion production. A more detailed report 
will appear elsewhere fll2[| . 
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FIG. 6. Order parameter \<f>\ versus et. 



VI. CONCLUSIONS AND OUTLOOK 

We considered the real-time dynamics of a coupled sys- 
tem of classical Bose fields and a quantized fermion field, 
represented by generalized mode functions, on a lattice 
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